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Abstract 

The density of accessible levels in the (^He,a7) reaction has been extracted 
for the ^^■^Dy, ^^^Er and ^''^Yb nuclei. The nuclear temperature is measured 
as a function of excitation energy in the region of - 6 MeV. The temperature 
curves reveal structures indicating new degrees of freedom. The heat capac- 
ity of the nuclear system is discussed within the framework of a canonical 
ensemble. 
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A challenging goal in nuclear physics is to trace thermodynamical quantities as functions 
of excitation energy. These quantities depend on statistical properties in the nuclear many 
body system and may reveal phase transitions. Unfortunately, it is difficult to investigate 
these aspects - both experimentally and theoretically. 

The density of levels as a function of excitation energy is the starting point to extract 
quantities like entropy, temperature and heat capacity. In the pioneer work of Bethe the 
level density was described within the Fermi gas model using a partition function for the 
grand-canonical ensemble. This picture of the nucleus as a gas of non-interacting fermions 
confined to the nuclear volume has later been modified. The phenomenological back-shifted 
Fermi gas model is a popular extension of the model, which simulates shell and pair 
correlation effects. This model works well at excitation energies where the level density is 
high, typically above the neutron binding energy. 

In addition to reveal statistical properties of nuclear matter, knowledge of the level den- 
sity is important in nuclear astrophysics. The level density is essential for the understanding 
of the nucleosynthesis in stars, where thousands of cross-sections have to be included in the 
calculations 0]. In parallel with these applications new theoretical approaches are emerging. 
Recent microscopic model calculations 0] include pair correlations as well as shell effects. 
With the recent shell model Monte Carlo method 0-0 one is able to estimate level densities 
in heavy nuclei up to high excitation energies. 

The theoretical progress has so far not been followed by new experimental data. The 
straightforward way to determine level densities is by counting discrete levels. However, this 
technique is restricted to light nuclei and/or low excitation energies where the experimental 
resolution is high enough to resolve individual lines in the spectra. In rare earth nuclei, the 
estimates obtained by counting levels are only valid up to ~2 MeV of excitation energy. 
A very useful experimental quantity is the average level spacing observed in slow neutron 
resonance capture. From these spacings rather accurate level densities for a certain spin 
value can be determined at the neutron binding energy region. In addition, the level density 
can be extracted from the shape of continuum particle spectra. However, if pre-equilibrium 
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particle emission takes place, this procedure may be doubtful since high-energy particles 
reveal high temperatures and, thus, too low level densities. 

Recently, the Oslo group has presented a new way of extracting level densities from 
measured 7-ray spectra PJ^. The main advantage of this method is that the nuclear system 
is very likely thermalized prior to the 7-ray emission. In addition, the method allows the 
simultaneous extraction of level density and the 7-strength function over a wide energy 
region. In this letter we report for the first time on experimentally deduced temperatures 
and heat capacities of rare earth nuclei in the 0-6 MeV excitation energy region. 

The experiments were carried out with 45 MeV '^He-projectiles at the Oslo Cyclotron 
Laboratory (OCL). The experimental data are obtained with the CACTUS multidetector 
array ||10[ using the (^He,a7) reaction on ^^^Dy, ^^"^Er and ^'^^Yb self-supporting targets. 



The charged ejectiles were detected with eight particle telescopes placed at an angle of 45° 
relative to the beam direction. An array of 28 Nal 7-ray detectors with a total efficiency of 
~15% surrounded the target and particle detectors. 

The assumptions behind the method of data analysis and techniques are described in 
Refs. and only a few comments pertinent to the present work are made here. The 

experimental level density is deduced from 7-ray spectra recorded at a number of initial 
excitation energies E, determined by the measured a-energy. These data are the basis for 
making the first-generation (or primary) 7-ray matrix, which is factorized according to the 
Brink- Axel hypothesis |[TT| , p!2| as 

PiE,E,) = aiE,)piE~E,). (1) 

From this expression the 7-ray energy dependent function a as well as the level density p is 
deduced by an iteration procedure [9], using the same Oth-order trial functions as described 
in Ref. [9], and the average slopes and the absolute density depend on this starting point. 
In the following we concentrate only on the level density and its fine structure, which is 
assumed to be independent of particular 7-ray decay routes. 

The extracted level densities for the ^^^Dy, ^^^Er and ^''^Yb nuclei are shown as data 



points in Fig. 1. The data for ^^^Dy and ^ Yb deviate slightly from the results previously 
published 0. In the present work the extraction procedure has been refined by omitting data 
with < 1 MeV, where the first-generation 7-spectra exhibit methodical uncertainties. By 
excluding these data points, the error bars for the resulting level density are reduced. 

The level density p{E) is proportional to the number of states accessible to the nuclear 
system at excitation energy E. Thus, the entropy in the microcanonical ensemble is given 
by 

S{E) = So + \np{E), (2) 

where, for convenience, the Boltzmann constant is set to unity {ks =1). The normalization 
constant 5*0 is not important in this discussion, since it vanishes in the evaluation of the 
temperature 

Small statistical deviations in the entropy S may give rise to large contributions in the 
temperature T. In order to reduce this sensitivity, the differentiation of 5* is performed by 
a least square fit of a straight line to five adjacent data points at a time. The slope of the 
straight line is taken as the differential of S at that energy. Thus, an effective smoothing 
of about 0.5 MeV is performed through this procedure. Since the energy particle resolution 
is around 0.3 MeV, this differentiation procedure will not significantly reduce the potential 
experimental information. 

The temperatures deduced are shown as discrete data points in Fig. 2. The data reveal 
several broad structures in the 1-5 MeV region, which are not explained in a Fermi gas 
description. The most pronounced bumps are located at 1.8 MeV and 3.2 MeV in ^^^Dy, 
at 2.5 MeV and 3.7 MeV in ^^^Er and at 1.8 MeV and 3.0 MeV in ^'^^yb. These structures 
are interpreted as being the breaking of nucleon pairs and, at higher energies, the possible 
quenching of pair correlations. 

The extraction of specific heat capacity is given by 
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(4) 



{dT/dE)v' 



which effectively means that the entropy S{E) has to be differentiated twice. From the 
scattered data points of Fig. 2, one immediately sees that this will not give meaningful 
information for Cy with the present experimental statistics. 

At this point one could introduce a strong smoothing of the temperature curve in order 
to extract Cy- However, it seems rather incidental how this should be done. One way to 
proceed is to introduce the canonical ensemble in the description of the nuclear system. Since 
the temperature enters as a fixed parameter in this formalism, the statistical uncertainty of 
T does not introduce any additional fluctuations. The cost of this nice feature is that one 
has to make an average over a wide excitation energy region. 

The partition function in the canonical ensemble 



is determined by the multiplicity of states at a certain energy. Experimentally, the mul- 
tiplicity corresponds to the level density of accessible states, p{En), in the present nuclear 
reaction at energy bin En- The widths of the energy bins are 120, 120 and 96 keV for the 
^^^Dy, i^^Er and ^'^^Yb nuclei, respectively. 

The mathematical justiflcation of Eq. (5) is that the sum is performed from zero to 
inflnity. The experimental information on p covers the excitation region 0-7 MeV, only. 
For the three nuclei investigated, the proton and neutron binding energies are ~8 MeV. For 
excitation energies above this binding energy it is reasonable to assume Fermi gas properties, 
since single particles may be excited into a region with a very high level density (continuum). 
Therefore, due to lack of experimental data above 7 MeV, the level density is extrapolated 
from 7 MeV to higher energies by the Fermi gas model expression [p]2[ 



where C is a normalization factor and a is the level density parameter. The best fits of ppc 
to data in the 3.5 - 6.0 MeV excitation region give a level density parameter of a =17.8, 
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(5) 



n=0 




(6) 
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19.0, and 18.7 MeV~^ for the ^^^Dy, ^^^Er and ^'^^Y nuclei, respectively. These values are in 
agreement with data from slow neutron resonances and the semiempirical formula for a 
with values in between A/IQ and A/S MeV^^. The theoretical level density functions are 
displayed as solid lines in Fig. 1 and are drawn from 3.5 MeV in order to visualize the fits 
to data. 

In order to determine how far the sum of Eq. (5) has to be performed, the probability 
density function 

P*^* = I pUE)e--^ITiE 
is shown in Fig. 3 for three typical temperatures. At low temperatures, say T = 0.3 MeV, 
the nucleons are scattered no higher than E =5 MeV. However, this upper limit increases 
very rapidly for higher temperatures. For the typical temperatures studied in this letter, E 
=30 MeV has been tested and found to be a sufficiently high upper limit for the summing 
in Eq. (5). 

The excitation energy in the canonical ensemble is given by the thermal average 

oo 

< E{T) >= Y: E^p{E^)e-''-'^ , (8) 

n=0 

where p is the level density (shown in Fig. 1) composed of an experimental and a theoretical 
part for the excitation regions below and above ~7 MeV, respectively. 

The smoothing effect implied by the canonical ensemble can be investigated by calculat- 
ing the standard deviation for the thermal average of the energy 



V< E"^ > - < E >2, (9) 



giving e.g. (Tb = 3 MeV at E =7 MeV. This again shows that the energy in the canonical 
ensemble is strongly smoothed for a given temperature. Thus, one cannot expect to find 
abrupt changes in the thermodynamical quantities. 

From Eq. (8) the temperature in the canonical ensemble can be studied as a function of 
< E >. In Fig. 2 the canonical ensemble gives a temperature dependence (solid lines) that 
coincides well with the average values found in the microcanonical ensemble (data points). 
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This gratifying behaviour of the canonical temperature encourages us to use the canonical 
ensemble to estimate the heat capacity Cy as well. The heat capacity can be deduced by 
simply calculating the increase in the thermal average of the energy < E > with respect to 
T 

CyiT) ^ (10) 

The deduced heat capacities for the ^^^Dy, ^^^Er and ^^^Yb nuclei as functions of < 
E{T) > are shown in Fig. 4. All nuclei display similar dependencies, reflecting that the 
thermal averages of the excitation energies smear out structures seen in the experimental 
level densities. The heat capacities are mainly following the theoretical values obtained by 
using Eq. (6) with the proper values of a deduced from the flts of Fig. 1. For comparison, 
we also show the simplified expression Cy = 1\faE assuming p ~ exp(2\/ aE) with a = 18.5 
MeV~^. The canonical heat capacity shows no traces of the fine structures found in the 
level density below ~4 MeV of excitation energies, as also was the case for the canonical 
temperature. 

In conclusion, for the flrst time temperature and heat capacity based on 7-ray spectra 
have been extracted for rare earth nuclei. The extracted temperature curves deduced for 
the microcanonical ensemble reveal structures of unknown origin. The structures in the 1 - 
5 MeV excitation region are interpreted as the breaking of nucleon pairs and quenching of 
the pair correlations. The heat capacity could only be extracted by the use of the canon- 
ical ensemble. These semiexperimcntal values, which are not expected to show any fine 
structures, are in agreement with the Fermi gas predictions. It would be very interesting to 
see if more realistic theoretical calculations can describe the fine structures observed in the 
microcanonical ensemble. 

The authors are grateful to E.A. Olsen and J. Wikne for providing the excellent ex- 
perimental conditions. We wish to acknowledge the support from the Norwegian Research 
Council (NFR). 
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FIGURES 
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FIG. 1. Extracted level density (points) for ^^^Dy, ^^^Er and ^"^^Yb. The error bars show the 
statistical uncertainties. The solid lines are extrapolations based on the Fermi gas model. The 

curves are in arbitrary units and separated with a factor of ~10 for better visualization. 
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FIG. 2. Observed temperatures as functions of the excitation energy E (data points with 
statistical error bars). The sohd hues are temperatures as functions of average excitation energies 
< E > deduced within the canonical ensemble. 
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FIG. 3. The Fermi gas model energy distribution p{E) of the canonical ensemble. 
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FIG. 4. The heat capacities Cy extracted within the canonical ensemble. The dashed curve 
displays the simplified Fermi gas expression for a =18.5 MeV~^. 
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